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1.  Introduction 

An  equivalent  aperture  excitation  method  (referred  to  as  B0R3  in  this 
report)  for  predicting  the  fields  coupled  through  a rotationally  symnetric 
aperture  in  a body  of  revolution  (BOR)  is  described  here.*  B0R3  has  an 
advantage  over  a conventional  near-field  scattering  BOR  code  [2]  (referred 
to  as  B0R2  in  this  report)  in  that  B0R3  is  characterized  by  greater 
sensitivity;  tiny  field  strengths  are  determined  more  accurately.  In  B0R2, 

, the  cavity  field  is  found  by  adding  the  impressed  and  scattered  field.  For 

small  field  levels,  such  as  those  associated  with  tiny  apertures  or  those 
at  locations  deep  within  a cavity,  this  computation  involves  a subtraction 
of  almost  equal  numbers  with  a corresponding  loss  in  accuracy.  In  B0R3,  this 
difficulty  is  avoided  because  the  cavity  field  is  computed  as  a radiation 
problem  once  an  equivalent  aperture  excitation  is  determined.  The 
improvement  in  accuracy  of  B0R3  over  B0R2  is  demonstrated  with  the  open- 
ended  cylinder  problem  discussed  in  Section  4. 

Based  on  the  method  of  moments  BOR3  is  an  extension  of  the  widely  used 
BOR  work  originating  at  Syracuse  University  [1,3]  . (A  computer  program 
implementation  of  B0R3  [9]  also  employs  a near-field  subroutine  developed 
at  Lawrence  Livermore  Laboratory  [2].)  The  details  of  the  Syracuse  BOR 
method  are  summarized  in  Appendix  A. 

There  are  two  principal  reasons  for  restricting  a study  to  BORs. 

First,  many  structures  of  interest  can  be  approximated  by  BORs.  These 
include  streamlined  bodies  such  as  aircraft  fuselages  and  missiles. 

* This  method  can  also  be  applied  to  non-rotationally-symmetric-aperture 
f BOR  problems  although  the  details  become  considerably  more  complex.  The 

general  method  is  described  in  [5]. 
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Second,  the  electromagnetic  characteristics  of  a BOR  conveniently  decouple. 
That  is,  a single  Fourier  mode,  or  circumferentially  varying  sinusoid,  of 
excitation  (impressed  field,  aperture  field,  etc.)  excites  only  the 
corresponding  mode  of  response  (far  or  near  field,  current  density,  etc.). 
This  permits  significant  savings  in  computer  time  and  memory.  With  modal 
decoupling  a matrix  formulation  of  the  problem  results  in  a number  of  small 
matrices  to  deal  with  rather  than  one  huge  matrix.  The  latter 
often  arises  with  use  of  an  arbitrary  surface  formulation. 

In  Section  2,  the  theory  behind  B0R3  is  described.  In  Section  3,  B0R3 
is  detailed  as  an  extension  of  the  Syracuse  BOR  method.  In  Section  4, 

B0R3  is  evaluated  by  comparison  with  B0R2,  classical  methods  and  experiment. 
These  methods  are  applied  to  both  spherical  and  finite-cylindrical  cavities 
with  rotationally  symmetric  apertures.  Finally,  in  Section  5 B0R3  is 
applied  to  a thin  circumferential  slot  in  a missile-shaped  body. 

A user's  description  of  a B0R3  computer  code  and  a Fortran  source 
listing  with  example  printout  are  available  [9] . This  version  of  B0R3  is 
limited  to  0-polarized  plane-wave  excitation. 

2.  Theory 

Figure  1 shows  a plane  wave  incident  upon  a perfectly  conducting  BOR 
with  a rotationally  symmetric  aperture.  The  surface  of  a BOR  may  be  formed 
by  constructing  any  contour  in  the  y-z  plane,  called  the  generating  curve, 
and  rotating  it  alaout  the  z-axis.  The  cylinder  and  conesphere  are  typical 
BORs.  The  basic  treatment  of  BORs  by  the  method  of  moments  has  been  set 


forth  by  Harrington  and  Mautz  [3]  for  radiation  and  scattering  problems. 

They  developed  a computer  code  for  this  method  which  has  proven  useful 
in  solving  a large  class  of  problems  [1,3],  Essentially,  the  task  is 
reduced  to  inverting  a set  of  matrices  each  of  order  equal  to  the  number 
of  expansion  functions  applied  to  the  BOR  generating  curve.  Each  matrix 
corresponds  to  a different  Fourier  component,  or  mode,  n,  of  cj.j.cumferential 
variation  of  fields  and  current  densities. 

The  basis  for  the  adaptation  of  the  Harrington-Mautz  (Syracuse)  BOR 
method  to  aperture  coupling  via  the  equivalent  aperture  excitation  method, 
B0R3,  is  an  equivalence  theorem  due  to  Schelkunoff  [4]  which  states  that  the 
original  problem  containing  sources  external  to  the  body,  can  be  replaced  with 
an  equivalent  problem  having  only  aperture  current  sources  for  excitation. 
Figure  2 illustrates  this  where  the  original  problem  depicted  in  (a)  is 
equivalent  to  the  sum  of  the  two  problems  in  (b)  and  (c) . In  (b) , the 
induced  surface  current  density  J°in  the  vicinity  of  the  aperture  is  found 
for  the  apertureless  (replace  aperture  with  conductor)  body  in  the 
presence  of  the  external  sources.  Then,  as  shown  in  (c) , = -J°  becomes 

an  apertiire  excitation  resulting  in  the  internal  fields  of  the  original 
problem.  For  BORs  both  problems  (b)  and  (c)  can  be  solved  by  proper 
modification  of  the  Syracuse  Method.  This  modified  method  is  B0R3.  It  is 
shown  in  [5]  that  BOR3  is  also  applicable  to  non-rotationally-symmetric- 
aperture  BOR  problems,  with  all  the  computational  advantages  of  BOR 
modeling  maintained,  although  the  details  are  then  considerably  more  complex. 

The  method  for  solving  the  problems  of  (b)  and  (c)  (Figure  2)  is 


outlined  as  follows: 


The  problem  depicted  by  Figure  2b  is  a straight-forward  BOR  scattering 
problem.  The  problem  of  Figure  2c  is  solved  by  first  replacing  the  current 
excitation  J with  an  equivalent  E-field  excitation.  The  latter  is  found 
by  dividing  the  aperture  into  thin  ribbon-like  cells,  or  ports,  along  the 
generating  curve,  short-circuiting  them,  and  exciting  each  one  in  turn 
with  a bi-directed  (since  the  surface  is  two-dimensional)  E-field,  each 
component  of  which  is  unity  across  the  excited  port  and  zero  elsewhere. 

The  circumferential  variation  of  each  port-E-field  excitation  is,  of  course, 
the  Fourier  component  or  mode  under  consideration.  Thus,  for  each  port 
E-field,  the  bi-directed  current  density  induced  at  every  port  is  easily 

obtained  by  solving  a straight-forward  BOR  radiation  problem.  This 

A A 

results  in  an  "aperture  admittance  matrix",  [y  ].  The  inversion  of  [y  ] 

and  post-multiplication  by  the  column  matrix  representation  of  J results 

in  the  desired  E-field  excitation.  The  latter  approximates  the  aperture 

tangential  E-field,  E , of  the  original  problem.  Excitation  of  the 

aperture  with  E is,  again,  a straight-forward  BOR  radiation  problem. 


3 . Computational  Details 

The  B0R3  method  is  an  extension  of  the  Syracuse  BOR  method.  The 
latter  is  detailed  in  references  [1,2,3].  A brief  summary  of  the 
formulations  and  notation  in  these  references  is  given  in  Appendix  A. 
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3.1  Aperture  Modelling 

The  aijcrture  is  considered  as  divided  into  P rotationally 
symmetric  ports.  For  simplicity  the  choice  of  location  and  extent  of  a 
port  is  constrained  to  contain  only  one  triangle  expansion  function  peak. 
The  triangle  functions  referred  to  are  those  which,  divided  by  the  BOR 
radius,  p,  are  cliosen  to  give  a piece-wise  linear  approximation  to  the 

current  density  along  the  generating  curve  of  the  apcrtureless  BOR. 

th  A 

The  n mode  aperture  E-field  excitation,  E , is  determined 

numerically  by  subsection  expansion  where  each  subsection  corresponds  to 

an  aperture  port.  Thus 


t = I 


p=l 


(1) 


where 


E = E e^"-^  X Nt);  + E ^ e^"^  X ‘'’(t)  I 
-p  p p p P ^ 


(2) 


Here  E ^ and  E are  the  n^^  mode  t and  d)  expansion  function  coefficients 
P P 

respectively  of  the  p^^  port  E-field,  and  X ^ (t),  X^*^  (t)  are  the  functions 

. . t 

which  describe  the  t dependence  of  these  expansion  functions.  Once  E^ 

and  E are  known,  the  cavity  field  con  be  found  by  solving  an  aperture 
P 

radiation  problem  as  outlined  in  Appendix  A.  The  remainder  of  this 

t d) 

section  deals  with  the  determination  of  E and  E 

P P 
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Note  that  a subscript  (or  superscript)  n is  omitted  from  the  symbols  of 


various  quantities  arising  in  this  development.  The  dependence  on  mode 
number,  however,  is  assumed.  Due  to  the  rotational  symmetry  the  derived 
quantities  for  each  mode  are  determined  independently.  Thus  this 
development  can  be  repeated  for  each  mode.  The  total  response  is  then  the 
sum  of  all  the  individual  modal  (positive  and  negative  n)  responses 
each  having  a (|)  dependence  given  by  the  complex  exponentials  in  (2) . The 
addition  of  a positive  and  negative  n response  exhibits  the  equivalent 
sin  (nij))  or  cos  (ncp)  dependence. 


3.2  Aperture  Admittance  Matrix 

The  elements  of  [y^] are  derived  here.  The  matrix  [y^J  relates 

A A 

column  vector  representations  of  E and  J [equation  (12)]. 


The  pq^*^  element  of  [y^ J is  determined  by  exciting  the  short- 
circuited  aperture  at  the  p^^  port  with  an  arbitrary  value  of  and 

recording  the  current  density  flowing  through  the  q port.  With  the 

* 

BOR  excited  by  the  generalized  voltage  excitation  column  vector  can 

be  partitioned  into  a t-component  sub-vector,  V , and  a <Ji-component 

P 


sub-vector,  V 


til  t 

The  j element  of  v has  the  form 

P 


V.^  = <w'^,E  > = 2itE 

J no  P p 0,p 


(3) 


Appendix  h reviews  the  formulations  (and  defines  the  symbolism)  for 
determining  a BOP,  current  density  given  an  impressed  or  aperture 
E-field. 
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similarly  for  the  (}i-component , 


V.*^=  <W  t , E > = 2ttE 
D n]  ' p p j,p 


(4) 


In  (3)  and  (4)  , 


F.^  = / pf.(t)  X ^(t)dt 

D»P  ■'gc  3 P 


(5) 


3»P 


= / 


pf  (t)  X ^ (t)dt 
gc  3 p 


(6) 


rr 


pf j (t)  = the  triangle  function 

the  weighting  functions  , and  tlie  symmetric  product  denoted  by 

"<  >"  are  defined  in  Appendix  A,  and  gc  stands  for  "generating  curve". 

Let  denote  the  triangle  function  with  peak  at  the  p^^ 

aperture  port.  Then  because  of  the  aforementioned  constraints  on  the 

extent  of  a port  the  functions  X ^(t)  and  X '*'(t)  are  zero  beyond  the  peaks 

P P 

of  the  j -1  and  j +1  triangle  functions.  Theref  ore , F . ^ and  F.*^  are  non-zero 
P P 3.P  3'P 

only  for  j = j^-  1,  j = j^,  and  j = jp^^'  (Further  considerations 
concerning  choice  of  xNt)  and  X^Ct)  are  given  later.)  Thus  the  excitation 
sub-vectors  have  the  form 


P 3p-l,F 

2 Ti  E ^F 

P 3p,P 

2 IT  E ^F^ 

P 3p+’'P 


V 


Le  "^Ft  - 

P 3p-l,P 

2TTE  *^Ft 
P 3p,P 

2rE  '^Ft  , 

P 3 +1»P 
0 ^ 


(7) 


(i 


t 

( 

i 


i 


c 


0 


Of  course  for  ports  located  at  the  junctions  of  generating  curve  and  axis 


(poles)  F.  -is  absent  for  j = 1 and  F.  , is  absent  for  j = nm  where 
Dp-1  ■’p  Dp+  1 -’P 

nm=  the  number  of  triangle  functions  on  the  generating  curve. 

The  "current"  flowing  through  the  port  due  to  a localized 

excitation  at  the  p^'^  port  can  be  obtained  from  the  generalized  network 

admittance  matrix.  The  latter  relates  and  ^ to  the  t and  (J>  generalized 

network  current  subvectors  and  respectively,  (Appendix  A)  by 


_ _p_ 

[Y<t>t]  ! ^ <P 


The  minus  sign  in  (8)  arises  since  V and  correspond  to  aperture 

P P 

(total)  E-fieids  and  not  "impiessed"  fields. 

For  simplicity,  the  correspondence  between  and  the  corresponding 

-*■  t <ti 

generalized  current  sub-vectors  is  made  only  at  t = t^,  the  location 

of  the  peak  of  the  triangle  expansion  function  pf.  (t) . Thus,  at  t = t^. 


; . J = ( i ) e^"*  = e^''^ 

'‘I  Pq  ^q 


. J = e^"^ 

-<5  Pq  <1 


where  I is  the  j element  of  I and  p is  the  BOR  radius  to  the 

3q  q q q 

point  t . 


y 


The  combination  of  ( 7,8,9)  and  removal  of  all  but  the 
rows  in  the  admittance  matrix  results  in 


tt 

t4> 

q»p 

V 

q»p 

4>t 

q*p 

^q»p 

. -2.  Y**’  . , , * y“*>  . F*?  * , f'"  , 

IT  ^ V ^q'=p  ^q'V  ^ V 


where  the  row,  6^^  column  element  of  [y^^l  and  ab=tt,  t(j),  (Jit,  and 

Ur  3 


The  aperture  admittance  relationship  for  all  P ports  then  becomes, in 
matrix  form. 


i = [y^J* 


where 


j = -i_ 


a = t or  and  the  q,p^^  element  of  [y^^] , ab  = tt,  t(|),  ijit,  or  is  given 
by  (11) . Note  that  j and  are  defined  as  in  (9)  but  with  J now 

q q -q 

representing  the  total  current  density  flowing  through  port  q (entire 

aperture  excited  as  opposed  to  only  E present) . 

~P 


3.3  Aperture  E-Field 

A 

The  aperture  tangential  E-field,  E , is  represented  by  (1,  2).  The 

coefficients  in  (2)  are  given  by  the  solution  to  (12)  for  e when  j is 

the  column  vector  representation  of  J . The  latter  is  known  from  the 

solution  to  the  aperturcloss  body  problem. 

The  X^(t),  X*^(t)  functions  are  chosen  to  be  essentially  pulses 
P P 

with  special  consideration  given  to  their  form  at  tfie  aperture  edges. 

One  pulse  is  associated  with  each  triangle  function  that  has  a peak 

within  the  aperture.  These  aperture  pulses  are  specified  by  the  corresponding 

triangle  expansion  function  numbers,  i,  and  the  lower  At^,  and  upper  AtT, 


1 
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extents  of  the  pulse.  Note  Figure  3.  The  constraints  on  the  pulse  width 
are  that 


At . < t . - t . , 

1 1 1-1 


At.  < t.^  - t. 

1 1+1  1 


The  <}i- component  of  E should  be  zero  along  a (j>  = constant  aperture 

edge,  and  the  t-component  becomes  infinite  there.  Therefore,  the  are 

P 

modified  by  setting 


At  ->■  At  /2  at  a lower  edge 

+ + 

At  ->■  At  /2  at  an  upper  edge 


and  the  are  modified  by  doubling  the  amplitude  of  the  portion  of  a 
pulse  closest  to  an  edge  (the  At  portion  of  a pulse  located  at  a lower 
edge  and  the  At^  portion  of  a pulse  located  at  an  upper  edge) . Thus  the 
t-coraponent  functions  at  aperture  edges  arc  really  step  pulses. 

Figure  3 exemplifies  these  pulses. 


4.  Evaluation  of  B0R3 

Cavities  of  convenient  shape  were  chosen  in  order  to  compare  B0R3 


with  classical  techniques.  Results  with  spherical  and  cylindrical  cavities 
are  presented  in  this  section.  Also  included  are  results  obtained  from 


both  the  conventional  application  of  a BOR  code,  B0R2,  and  experiment. 
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Figuz'C  3.  IVplcal  (a)  Triangle  Functions  for  Current  Density  Expansion 
Along  BOR  Generating  Curve  and  Modified  Pulse  Functions  for 
Aperture  E-field  (b)  t-directed  Component  X and  (c) 

♦“directed  Component  X^.  ^ 

P j 

j 

I 


The  purpose  of  these  comparisons  was  twofold.  First  was  the  determination 
of  the  number  of  triangle  functions  per  aperture,  "nam",  required  for 
reasonable  accuracy  in  field  (~20%).  Second  was  the  evaluation  of  the 
r effect  of  severe  size  variation  in  adjacent  triangle  segments  (four 

^ segments  per  triangle  extent  [1]).  These  considerations  are  pertinent 

if  the  height  of  the  circumferential  aperture  is  only  a tiny  fraction  of 
the  cavity  length.  For  a given  nam  and  a reasonable  number  of  triangle 
^ functions  on  the  cavity  generating  curve,  "nm",  a large  variation,  e.g.  10:1, 

in  adjacent  segment  size  may  result. 

f 4.1.  Spherical  Cavity 

The  equivalent  aperture  excitation  method,  B0R3,  was  applied  to  the 
spherical  cavity  aperture  coupling  problem  studied  by  Senior  and 
Desjardins  [6].  The  latter  employed  a spherical  mode  matching  method. 

The  case  of  a ka  = 2.5  sphere  with  30°  half-angle,  front-end  aperture  and 
excited  by  a unit  amplitude,  linearly  polarized.  E-field  incident  along 
the  axis  and  toward  the  aperture  side  of  the  body  is  shown  in  Figure  4. 

The  magnitude  of  the  internal  E-field  along  the  axis  is  plotted. 

Only  the  n = ±1  mode  (sin(J,  cos4>)  was  excited  because  the  impressed  plane 
wave  was  axially  incident.  The  BOR3  result  for  a total  of  nm  = 15  triangle 
functions  along  the  generating  curve  of  the  apertureless  sphere  is  in  good 
agreement  with  the  spherical  mode  matching  result.  Five  of  these 
triangles  traversed  the  aperture  (nam  = 5).  This  necessitated  a 


SPHERICAL  MODE  MATCHING 
(SENIOR  AND  DESJARDINS) 


(BACK)  (CENTER)  (APERTURE) 


three-to-one  variation  in  segment  size  with  no  appreciable  loss  in 
accuracy.  The  restriction  of  uniform  segment  size  and  nm  = 11  (not  plotted) 
resulted  in  only  one  triangle  peak  lying  in  the  aperture.  However, 
the  aperture  field  then  differed  from  the  nm  = 15  case  by  about  6 db. 

The  result  with  nm  = 13,  nam  = 3 (Figure  4)  differs  by  less  than  20%. 

Thus  a minimum  of  three  triangle  functions  should  "excite"  the  aperture. 

It  is  not  surprising  that  one  aperture  expansion  function  is  inadequate, 
even  for  tiny  aperture  heights.  Small-Hole  theory  predicts  that  two  modes-  a 
"magnetic"  (proportional  to  J^)  and  an  "electric"  (proportional  to  V • J^)  - 
are  generally  significant. 

In  Figure  5 the  field  along  the  axis  of  a ka  = 1.0  sphere  with  10° 
half  angle  aperture  is  plotted.  The  results  for  a 10:1  segment  size 
variation  are  in  reasonable  agreement  with  those  for  a 2:1  segment  size 
variation.  In  Figure  6 the  magnitude  of  the  scattering  current  densities 
on  the  corresponding  apertureless  sphere  are  plotted  for  both  cases.  The 
most  significant  discrepancy  is  in  the  (J)-component  of  current  density. 

Also  plotted  in  Figure  5 is  the  spherical  mode  matching  result.  Because 
of  the  smallness  of  the  sphere  and  tinyness  of  the  aperture  it  is  difficult 
to  determine  which  method  is  the  more  accurate;  the  spherical  mode  matching 
method  suffers  from  slow  convergence  especially  at  the  aperture  edges. 
Apparently  a 10:1  variation  in  adjacent  segment  size  is  not  very  reliable 
whereas  a 3:1  variation  is. 
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SPHERICAL  MODE  MATCHING 


Figure  5.  Amplitude  of  E-Field  Along  Axis  of  Sphere,  ka 
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4.2  Cylindrical  Cavity 

Further  evaluation  of  B0R3  was  obtained  with  application  to 
an  open-ended  cylinder  problem.  The  comparison  was  with  B0R2,  classical 
theory,  and  experiment.  The  magnitude  of  the  E-field  along  the  axis  of  the 
cylinder  is  plotted  in  Figure  7.  The  excitation  v/as  a linearly  polarized 
unit  plane  wave  axially  incident  on  the  opci:-ended  side.  Again,  only 
the  n = + 1 mode  was  excited  due  to  the  axial  incidence.  The  conventional 
application  of  the  BOR  code,  B0R2,  and  the  equivalent  aperture  excitation 
method,  B0R3,  agree  well  to  a distance  of  about  . 2A  into  the  cylinder.  The 
scale  factor  difference  between  them  is  attributed  to  the  different  ways 
in  which  each  models  the  aperture;  an  increase  in  the  number  of  expansion 
functions  should  bring  both  methods  into  better  agreement.  Beyond  . 2X 
into  the  cylinder  the  B0R2  result  tends  to  level  off  and  approach  the 
"closed  body"  (apertureless)  result  whicli  is  also  plotted.  This  is  to 
be  expected  since  the  "closed  body"  field  is  essentially  the  error 
field  (ideally  zero)  resulting  from  the  approximations  in  method-of-moment 
modelling  and  the  errors  associated  with  the  aforementioned 
(Section  1)  addition  of  impressed  and  scattered  fields.  B0R3  is  not 
limited  by  the  "closed  body"  field.  Thus  it  demonstrates  a minimum  of 
perhaps  20dB  improvement  in  sensitivity  over  B0R2. 

A check  on  the  B0R3  result  was  obtained  from  a classical 
v/aveguide-mode  analysis.  The  modal  components  were  determined  by  requiring 
zero  tangential  E-field  on  the  cylinder  surface  except  in  the  aperture  where 
the  E-field  predicted  by  B0R3  was  prescribed.  Five  electric  and  five 
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Figure  7.  Amplitude  of  E-Field  Along  Axis  of  an  Open-Ended  Finite  Cylinder. 


magnetic  modes  were  used.  Due  to  the  thinness  of  the  cylinder  (the 
frequency  was  well  below  cutoff)  the  dominant  mode  alone  essentially 
resulted  in  the  fields  beyond  the  aperture.  As  shown  in  Figure  7,  there 
is  no  observable  difference  between  the  cylindrical  mode  matching  method 
and  B0K3 . 

A further  check  on  the  computational  methods  was  obtained  with 
an  anechoic  chamber  experiment.  An  electrically  short  diode-loaded 
probe  and  a high  resistance  transmission  line  were  used  to  minimize 
measurement  apparatus  interference.  The  probe  was  mounted  on  a non- 
metallic  "plunger"  with  placement  holes  radially  spaced.  The  plunger 
was  mounted  on  a calibrated  dielectric  rod.  In  this  v/ay,  all  positions 
of  the  probe  could  be  determined  accurately.  A sampling  of  the  results  are 
plotted  in  Figure  7.  They  are  adjusted  to  account  for  the  dielectric  support 
which,  experiment  showed,  increases  the  observable  field  by  about  20%. 

There  is  good  agreement  between  theory  and  experiment  for  distances  of 
about  .15X  or  less  into  the  cylinder.  For  greater  distances  the 
experimental  results  were  not  reliable,  and,  therefore,  they  are  omitted 
from  the  graph.  The  reason  for  this  is  that  at  the  . 15X  point  the 
dominant-mode  field  was  approximately  30dB  below  the  impressed  field,  and 
harmonics  of  the  fundamental,  which  were  closer  to  or  above  the  cut-off 
frequency  of  the  "guide,"  became  significant.  Attempts  at  filtering 
the  transmitter  to  suppress  these  harmonics  could  not  appreciably  extend 
the  curve  due  to  the  severe  attenuation  of  the  dominant  mode  at  the 
fundamental  frc!quency. 
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5.  Thin  Circumferential  Aperture 


The  equivalent  aperture  excitation  method,  B0R3,  was  applied  to  the 
missile-like  structure  shown  in  Figure  8.  The  aperture  is  modelled  simply 
as  a small  separation  between  two  sections.  It  is  noted,  however,  that 
actual  apertures  are  often  considerably  more  complex  since  the  section  walls 
can  overlap  at  the  junction  as  indicated  in  Figure  14,  or  the  aperture 
can  be  filled  with  a dielectric  material  or  sealed  with  a conductive  putty. 
An  extension  of  the  foregoing  technique  is  available  which  can  account  for 
the  aperture  complexity.  This  is  discussed  later  in  this  section  although 
computations  have  not  as  yet  been  performed. 


5.1  Oblique  Incidence  Plane  Wave  Excitation 

Figures  9-13  show  the  E-field  in  the  vicinity  of  the  aperture  for  a 

unit,  0-polarized  impressed  field  arriving  from  the  0 = 45°  direction  in  the 

x-z  plane.  Each  figure  plots  three  lines.  Each  line  corresponds  to  a 

path  indicated  by  one  of  the  dashed  lines  shown  in  Figure  8.  The  azimuth 

orientation  of  each  path  corresponds  to  that  of  peak  amplitude.  The  p- 

and  z-  components  of  E-field  (E*^,  E*')  each  vary  as  cosnij),  and  the  4>-component, 

p z 

E^,  varies  as  sinncj) . Plots  are  shown  for  the  n = 0 and  1 modes  except  for 
0 

E,  which  is  identically  zero  due  to  the  0-polarization  of  impressed  field. 

w 

Other  modes  are  not  presented  since  the  n = 2 modal  fields  were  found  to  be 

more  than  20db  below  the  n = 0,1  fields.  This  is  not  surprising  since 

the  number  of  significant  modes  is  limited  to  all  n such  that 

|n|  ^ 1 + C|^  where  is  the  largest  circumference  of  the  body  in  wavelengths 

[1].  However,  it  is  interesting  to  note  that  the  n = 1 mode  is  significant 

even  though  the  structure  is  quite  thin.  Although  the  sinusoidal 
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circumferential  variation  of  the  n = 1 mode  would  make  it  insignificant  for 


far-field  scattering  it  is  certainly  of  significance  in  aperture  coupling. 

As  expected,  exhibits  a dipole-like  behavior  (Figure  10) . For  p 
close  to  an  aperture  edge,  where  the  charge  density  is  infinite, 

P 

0 

becomes  large.  Also  the  phase  of  E^  changes  rapidly  in  passing  by  the 

0 1 

aperture.  The  behavior  of  E^  and  E^  are  also  as  expected  in  that  close  to 
the  cavity  surface  they  are  small  (except  at  the  aperture)  (Figs.  9,  11) 
and  close  to  the  axis  E^  is  small. ( From  symmetry,  E^  should  vanish  along 
the  axis. ) 

The  number  of  triangle  functions  used  in  computing  these  results  was 
nm  = 44.  (The  corresponding  generalized  impedance  matrix  was  of  order  88 
since  the  current  density  is  bi-directed. ) Also,  as  suggested  by  the 
results  of  Section  4,  three  triangle  peaks  were  chosen  to  lie  in  the  thin 
aperture.  This  resulted  in  a 3:1  triangle  segment  size  variation  near  the 
aperture. 


5.2  Complex  Aperture  as  Equivalent  Load 

In  the  beginning  of  this  section  it  was  pointed  out  that  in  practice 

apertures  are  usually  considerably  more  complex  than  simple  circumferential 

openings.  A typical  aperture  is  depicted  in  Figvure  14.  A method  for 

including  the  effect  of  the  aperture  shape  and  material  composition  in 

computing  cavity  fields  is  presented  here.  Essentially  it  involves  the 

determination  of  an  "aperture  load"  matrix,  [y  ] , which  can  be  combined 

x« 

A A 

with  [y  ] such  that  e,  the  column  vector  representation  of  E , can 
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Figure.  14.  Typical  Circumferential  Aperture. 
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readily  be  determined. 

As  discussed  in  Section  2,  both  j,  the  current  vector  representation 
A A 

of  J , and  also  [y  ] are  computed  with  the  aperture  "short  circuited". 

A 

Hence,  j and  [y  ] are  independent  of  the  aperture  complexity  (or  "loading"). 

A 

With  the  aperture  present  and  excited  with  J (Figure  15a)  the  problem  can 
be  represented  in  networ)c  form  as  shown  in  Figure  15b.  The  generalization 
of  (12)  to  include  the  aperture  complexity  is  then 

j = ^[y^]  + 

A A A 

Once  [y  ] is  )cnown,  e and  thus  E , can  be  determined  from  (14).  With  E 

)cnown  the  internal  fields  can  be  computed,  as  before,  from  a straightforward 

aperture  radiation  problem  (Figure  15c). 

For  thin  apertures  the  load  matrix  elements  relating  the  t-directed 

components  of  current  and  E-field  can  be  approximated  with  quasi-static 

solutions.  The  shape  of  the  conducting  shell  at  the  aperture  presents 

essentially  a capacitance  the  value  of  which  can  be  approximated  by 

Contour  Mapping.  Also  the  conductance  and  susceptance  of  material  (other 

than  free  space)  filling  the  aperture  can  be  determined  by  mapping  methods. 

The  ()i-directed  components  of  current  can  easily  s)cirt  the  thin  aperture. 

Thus  only  the  t-directed  components  of  loading  need  be  considered. 


y‘ 


f) 


(14) 
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6.  Conclusions 


A moment  method  for  predicting  aperture-coupled  fields  within  a 
rotationally  symmetric  cavity  has  been  presented.  This  method,  B0R3, 
employs  an  equivalent  aperture  excitation  concept  which  demonstrates 
greater  sensitivity  over  more  conventional  procedures  if  the  aperture  is 
thin.  This  improvement  is  important,  for  example,  in  analyzing  the  fields 
coupled  through  thin  circumferential  apertures  in  missile-like  structures. 
A computer  progrcim  with  user  manual  is  available. 
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APPENDIX  A 


REVIEW  OF  METHOD  OF  MOMENT  MODELLING 


OF  A BODY  OR  REVOLUTION 


The  necessary  relationships  for  the  Harrington-Mautz  [l,3]  method- 
of-moment  modeling  of  a body  of  revolution  (BOR)  are  briefly  described  here. 
Only  conducting  bodies  arc  considered  although  surface  formulations  for 
material  bodies  exist  to  which  this  theory  can  be  extended  [?].  This 
method  is  based  on  the  E-field  formulation.  Thus  both  radiation  and 
scattering  problems  are  solved  by  a single  procedure,  and  also  the  thin- 
walled  cavity  problem  can  be  handled.  The  latter  is  important  to  aperture 
coupling  considerations. 

The  E-field  formulation  relates  a free-space  surface  current  density 
J to  the  tangential  component  ot  E along  the  surface  S of  a conducting  body. 
In  the  case  of  scattering  by  a perfect  conductor 

e|  =E^L  +E®|.  =0 

~ ' tan  ' tan  ~ ' tan 

X s 

along  S where  E is  the  impressed  field  (body  absent),  E is  the 

field  radiated  by  J,  and  denotes  "component  tangential  to  S". 

Hence,  the  E-field  formulation  becomes 


5 Itan  = 


where 


j(4j//j  ds  + V~  (V‘J)e^ 


(or,  S 
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Here  R is  the  distance  between  source  and  field  points  on  S,  the  integrations 
are  v/ith  respect  to  the  source  points,  and  the  remaining  symbols 
have  their  usual  meanings.  In  the  case  of  radiation  is  known  along 

S and  the  E-field  formulation  then  becomes 

-e|  = L(J)  (A-2) 

An  approximation  to  J satisfying  ^A-1)  can  be  found  by  the  method 

of  moments.  Consider  the  coordinate  system  of  Figure  1.  A set  of  expansion 

functions  = f . (t)e^*^'^t  and  f'^  . = f . (t)  (t  and  $ are  unit  vectors) 

~nj  3 -n]  j 

and  the  syr metric  product 

< J,  , fl  • J2  ‘>‘>♦'5' 

are  chosen.  (The  "gc"  stands  for  BOR  generating  curve.)  The  f ^ (t) , which 
express  the  t variation  of  the  expansion  functions,  are  selected  in 
accordance  with  the  method  of  subsections  [s].  In  the  Harrington-Mautz  BOR 
model  they  are  triangles  (overlapping)  divided  by  the  BOR  radial  coordinate 
p.  This  choice  permits  differentiation  of  the  expansion  functions  and  a 
non-zero  value  where  the  poles  of  the  BOR  meet  the  axis.  The  Fourier 
components  are  chosen  so  that  the  modal  decoupling  property  of 

rotationally  symmetric  bodies  applies.  Thus  for  nm  triangles  on  the 
generating  curve  J is  expanded  as 

J = I I (I  ^ f . (t)e^”'^t  + l'’’.f . (t)e^"'^^)  (A-3) 

^ ^ n:  D nj  3 

n=-oo  jr=l 
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where  1 and  I vire  the  unknown  coefficients  to  be  determined  such  that 
nj  nj 


(A-1)  [or  (A-2)]  is  satisfied.  The  inner  product  of  both  sides  of  (A-1) 
arc  taken  with  respect  to  each  of  the  weighting  functions  = f^(t)e 

and  W t = f . (t)e  resulting  in  the  matrix  equation 


V = [z]l 

n n-*  n 

for  each  n.  Upon  partitioning,  (A-4)  becomes 

_tti  I r„t(ti  1 r 


1 

[ 

r 

< 

ft 

! 

[z^^]  1 

n 

n 

1 

n 

< 

3 -©- 

1 

1 

n 

VJith  a=t  or  ^ the  i^^  element  of  each  V ^ and  I ^ subvector  is 

n n 

< W ? , > and  I ? respectively.  The  i,j^^  element  of  each 

ni  ni  , '-n-' 

b 

submaLrix  is  < , L > where  ab  = tt,  (}it,  t(J> , (|)(j). 

Each  mode  of  J is  determined  independently  from  (A-3, 4) resulting  in 
significant  savings  in  computer  processing  time  and  memory.  The 
solution  to  (A-4)  can  be  expressed  as 


I = [y  ] V 
n n n 


where  [y  ] = [z  ] ^ is  the  generalized  admittance  matrix, 
n n 

The  n^^  mode  radiated,  or  scattered,  field  E can  be  determined 

-n 

from  (A-5)  by  reciprocity.  A current  element  II  is  located  at  the  field 
poir.t  cf  interest,  r,  and  oriented  parallel  to  the  polarization  of  interest. 
If  is  the  mode  excitation  vec;tor  due  to  11^  then 


*f  ^ a 

Here  is  the  transpose  of  V , V = , the  j element  of  V is 

n n n - - x n 

n 

W * 

cl  r r X 

< F .,E  > where  a = t or  <ii,  and  E is  the  field  radiated  by  11  . Both 

n: 

far  [if  3]  and  near  [2] fields  can  bo  determined  from  (A-6) . 


37 


RKFERl'i'JCES 


[ 

I 

t 

t 


I 


[1] .  R.F.  Harrington  and  J.R.  Maut^,  "Radiation  and  Scattering  from 

Bodies  of  Revolution,"  Final  Report,  prcjiarcd  for  Air  Force 
Cambridge  kesenrcii  Laboratories,  Bedford,  Mass.,  under  Contract 
No.  F1‘j6^8-67-C-0233,  duly  ]y69. 

[2] .  R.M.  Bevensee,  "S3F  - SYP/Ll.Ll  The  Syracuse  Computer  Code  for 

Radiation  and  Scattering  fiom  Bodies  of  Revolution,  Extended  for 
Near-Field  ComiJU tat  ions , " RejDort  No.  L'CRL-51622,  prepared  for 
U.S.  Atomic  Energy  Commission  under  Contiact  No.  W-7405-Eng-48 
with  the  Lawrence  Livermore  Laboratory,  University  of  California, 
Livermore,  CA,  May  1974. 

[3] .  J.R.  Mautz  and  R.F.  Harrington,  "Generalized  Network  Parameters  for 

Bodies  of  Revolution,"  Scientific  Report  No.  1,  Contract  No. 
F-19628-67-C-0233 , Air  Force  Cambridge  Research  Laboratories, 
Bedford,  Mass.,  with  Syracuse  University,  Syracuse,  N.Y.,  May  1968. 
Or  see  instead:  J.R.  Mautz  and  R.F.  Harrington,  "Radiation  and 

Scattering  from  Bodies  of  Revolution,"  Applied  Scientific  Research, 
Vol.  20,  1969. 


I 


[ 

I 


[4] .  S.A.  Sohelkunoff,  "Field  Equivalence  Theorems,"  Comm,  on  Pure  and 

Applied  Math.,  Vol.  4,  pn.  43-59,  June  1951. 

[5] .  H.K.  Schuman,  "Coupling  Through  Non-Azimuthally  Symmetric  Apertures 

in  a Body-of-Revolution" , Internal  Memorandum,  Dept,  of  Electrical 
and  Computer  Engineering,  Syracuse  University,  July  1974. 

[6] .  T.B.A.  Senior  and  G.A.  Desjardins,  "Electromagnetic  Field  Penetration 

into  a Spherical  Cavity,"  IEEE  Trans,  on  Electromagnetic 
Compatibility , Vol.  EMC-16,  pp.  205-208,  November  1974. 

[7] .  Private  Communications  with  1)L.L.  Tsai  and  T.K.  Wu  at  the  University  of 

Mississippi,  Oxford,  Mi.ssiru;  Ippi  , and  2)  R.F.  Harrington  and  J.R. 

Mautz  at  Syracuse  University,  Syracuse,  New  York. 

[9].  R.F.  Harrington,  Field  Ccinpur-^tion  by  Moment  Methods,  The  Macmillan 
Co.,  New  York,  1968. 


[9].  H.K.  Schuman,  "Coupling  tlirough  Rotational ly  Symmetric  Apertures  in 

Cavities  of  Revolution  - Computer  Code  (B0R3)  Description",  Tech  Memo, 
Electrical  and  C miuter  Engii.  cring  Dept.,  Syracuse  University, 
Syracuse,  NY  13-10,  September,  1‘,'76. 


38 


